Theory for wavelength-resolved photon emission statistics in single-molecule spectroscopy 
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We derive the moment generating function for photon emissions from a single molecule driven by laser 
excitation. The frequencies of the fluoresced photons are explicitly considered. Calculations are performed for 
the case of a two level dye molecule, showing that measured photon statistics will display a strong and non- 
intuitive dependence on detector bandwidth. Moreover, it is demonstrated that the anti-bunching phenomenon, 
associated with negative values of Mandel's Q-parameter, results from correlations between photons with well 
separated frequencies. 
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Single-molecule spectroscopy (SMS) Ol] provides a de- 
tailed glimpse into our natural world. Typically, SMS exper- 
iments rely upon broadband detection of fluoresced photons 
to monitor molecular dynamics, however in studies of reso- 
nance energy transfer (RET) l2] and semiconductor quantum 
dots iit] considerably more information can be obtained by re- 
solving photon emission into two color channels. Recently 
10], multi-channel detection schemes have been introduced to 
extend the capabilities of SMS still further 

The information obtained by SMS is useful only if it can 
be readily interpreted. SMS has received considerable theo- 
retical attention (see reviews jl] and references within), but 
most of this work ignores any consideration of photon color 
A notable exception is the treatment of RET, which has been 
considered in detail |@]. A related treatment of "frequency re- 
solved" photon counting, including quantum evolution of the 
molecule, has also been proposed by us 10] ■ However, these 
studies rely upon a direct correspondence between individ- 
ual spectral transitions and experimental detection channels. 
This picture may be adequate for well resolved transitions and 
certain experimental conditions, but falls short of providing a 
complete theoretical description of emission spectroscopy at 
the single-molecule level. 

Within the field of quantum optics, time correlations be- 
tween spectrally resolved photons have been studied both ex- 
perimentally 101 and theoretically for the case of resonance 
fluorescence from 2-level atoms. These studies also rely upon 
a direct correspondence between individual spectral transi- 
tions (in the dressed-atom picture |10]) and the frequency of 
the emitted photons to enable elementary interpretation of ex- 
periment and simplified theoretical analysis. This letter intro- 
duces a general formalism to describe single molecule photon 
emission that does not presume simplifying characteristics of 
the molecular system or detection apparatus. Our results may 
be directly applied to model systems and lay the groundwork 
for development of controlled approximation schemes in the 
study of more complex condensed-phase systems. 

In previous work, we (T^ and others [l2ll have intro- 
duced the generating function formalism for calculation of 



single-molecule photon counting statistics without spectral 
resolution. Such broadband photon statistics may be calcu- 
lated by monitoring the number of times that spontaneous 
emission occurs as the molecule evolves. Within the Marko- 
vian limit for molecular dynamics, spontaneous emission is a 
simple rate process and these emission events may be treated 
purely classically, even though the underlying dynamics may 
involve facets of quantum evolution. Calculation of photon 
counting moments proceeds via introduction of the generating 
function for spontaneous emission events 0(5, = (s"'^*'') 
where n{t) is the number of emissions in the interval [0, t] 
and the factorial moments of this quantity follow immediately 
by differentiating G with regard to the auxiliary variable s 
and evaluating at s = 1. The equations of motion for G{s, t) 
(and by extension the factorial moments) involve only mini- 
mal complications beyond the usual quantum master equation 
approach used to solve for density matrix dynamics lllll . 

In contrast to the above, if the frequency of emitted photons 
is measured, it becomes impossible to proceed via simple clas- 
sical arguments. Decay of an electronic excitation into a par- 
ticular field mode or narrow subset of modes can not be moni- 
tored by simply counting instantaneous spontaneous emission 
"events"; such a process is fundamentally non-Markovian. 
However, the definition of the generating function may be 
generalized to allow for calculation of factorial moments with 
frequency resolution by explicitly introducing a quantum me- 
chanical description of the radiation field. We take 



G{s,t) 



exp 



^ln(ske)aj,^(i)ake(i) 



ke 



(1) 



A/'cxp 



k6 



*Present Address: Center for Nonlinear Studies, Los Alamos National Labo- 
ratory, Los Alamos, NM 87545 USA 



Here, the averaging operation has its usual meaning (...) = 
Trace{...p(0)} involving the initial density matrix and a full 
trace over both the fluorophore and radiation field degrees 
of freedom. Creation and annihilation operators for pho- 
tons with wavevector k and polarization e have been in- 
troduced to express s"^*^ from the broad-band definition as 

exp[ln(s)X;ue^ke(i)] with N^,{t) = aL(*)ake(i) rep- 
resenting the Heisenberg picture number operator for each 
mode. The generalization from s to s has been made to facil- 
itate extraction of spectral information. The second equality. 
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involving the normal ordering operator M, follows from stan- 
dard operator identities ifisll . Taylor expanding both expres- 
sions around Ske = 1 reveals that the multivariate factorial 
moments of the number operators A^ke are obtainable by the 
traditional differentiation rule at Ske = 1 and that these mo- 
ments are most conveniently expressed as a normally ordered 
product of creation and annihilation operators for each mode 
appearing in a given moment. For example, we find 



5"+™G(s,i) 



'^*ke'^*k'e' 



(2) 



s=l 



([(«L)"(«Ls')"(ak's')™(«ke)"](t)> 



with the expected generalization applying to moments involv- 
ing more than two modes. The above introduces the notation: 
iV(™) = N{N - 1) . . . (A^ - m + 1). 

To make further progress, we specify the form of the Hamil- 
tonian governing the time evolution of the operators discussed 
above 111. 



H(t) = H,+Hr + Hi-[ D+{fio ■ El 



h.c. 



(3) 

Hs is the Hamiltonian for the system (atom or molecule) of 
interest, which will always be modeled with two electronic 
states (ground and excited) coupled to nuclear degrees of free- 
dom. Hfj = X^k- ^kflk^flkE is the Hamiltonian for the quan- 
tum radiation field {ujk ~ ck with c the speed of light). The 
last term in parentheses reflects a semi-classical coupling be- 
tween the applied laser field (assumed monochromatic with 
frequency lul and amplitude El) and the system within the 
dipole approximation (D = ij,q{D^ + D^) is the dipole mo- 



ment operator for the system consisting of terms that raise (+) 
and lower (— ) the electronic state of the system) and rotating 
wave approximation (RWA) IToll . Hj describes the interac- 
tion between the system and the modes of the quantized elec- 
tromagnetic field, also within the RWA and dipole approxima- 
tion 



iJ/ = 7^ (-i(£ • n^)D+a\^e + h.c.) 



(4) 



ke 



In the above 7 = y where ujq is the transition frequency 
between excited and ground electronic states 13, e is the per- 
mittivity and V is the volume of the cubic box used to quan- 
tize the field {V ^ oo below and does not appear in any final 
results). 

The Heisenberg equations of motion for the creation and 
annihilation operators evolving with dynamics dictated by eq. 
|3]may be formally integrated to yield lUOll 



e 



ak. (0) 



(e • iJia)D- (r) e^'^'-^^dT 



and the conjugate expression for a]^^{t) 



(5) 

For later conve- 
nience, we have introduced the slowly varying rotating frame 
operators D^{t) = D^{t)e^^'^'^* and have set cukL = ^^k — 
LOL- From this, it is readily seen that ake(i) commutes with 
ake(t) and similarly for a^^(t) and a\^^{t). This fact, along 
with the assumption that the initial time total (system and 
radiation field) density matrix is a direct product between 
the system and the vacuum state for the field ( i.e. p(0) = 
o-s(O) (8> |0) (0|) allows us to reformulate eq. [Das IH] 



t t 

G(t,V) = (r^exp|^^(ske-l)y j{e-iJi^fb+{u)D-{v)e-^^''-^:--"Uudv\ ). (6) 





The operator TV acts on all operators to the right of it by first 
arranging all "H-" operators to the left of all "-" operators and 
subsequently placing all "-" operators in standard time order 
(latest times at the left) and all "H-" operators in reversed time 
order (latest times at the right). The advantage of eq. |6]over 
either expression in eq. [T]is that the generating function is now 
defined solely in terms of the evolution of the system, which 
allows us to pursue actual calculations as detailed below. 

Eq. [U provides a theoretical route toward arbitrary photon 
counting moments. For simplicity and to make connection 
with possible experiments, we specialize to the case that pho- 
ton detection is insensitive to propagation direction and polar- 
ization of the emitted photons and also assume that the detec- 
tors have finite resolution, registering the arrival of all photons 
within a window of width A around a central frequency lo. We 



define a number operator for photons within this window 

Ni.., = V Nv. (7) 



ke:(w-A/2)<Wfe<(a)+A/2) 



Combining the above definition with eqs. |2]and|6]and pro- 
ceeding to the continuum limit {V oo) leads to the conclu- 
sion that 



To 

\ uj-A/2 



(8) 



duji D+ {u)D- {v)e-"^'^^''-''^dudv 



where To = 3° e^°J:i ■ Eq. [8] applied to the case m = 1 counts, 
on average, the number of photons within a given frequency 



window emitted in time t by the externally excited molecule. 
The time derivative of this quantity evaluated in the t ^ oo 
limit reproduces the usual expression lflol[l6ll for the spectrum 
of fluoresced radiation. Also, in the limit that A ^ oo, eq. [8] 
reduces to Mandel's expression ITtIi for the factorial moments 
of photon emission as detected in broadband measurements 
(i.e. no frequency resolution). 

The Tj^ operator in eq. [8] insures that all correlation func- 
tions appearing in (ivf'"-' 



' (w A) f^^)) ^® °^ '■^^ form 

. . . D+{Ura)D-{Vm) ■ ■ ■ D~ (vi)) 
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with U,n > 



-1 > . ■ .ui and v„i > > . . .vi. Cor- 



relation functions with such time ordering may be calculated 
within in the Markov limit for system dynamics ifioll via an 
extension of the quantum regression theorem ifisll . 

It follows that the explicit calculation of moments in eq. [8] 
is straightforward in principle, involving only diagonalization 
of the rotating-frame evolution operator for system dynamics 
and elementary integrals over time and frequency. In practice, 
however, the procedure is complicated and will be specified in 
detail elsewhere ifisll . 

For concreteness, we present predictions for the low tem- 
perature spectroscopy of a single 2-level dye molecule. The 2- 
level molecule is specified by i/s = {hujo/2) {\e){e\ — \g){g\) 
and = {D ~ l5)(e|) with e and g designating 

the excited and ground states. Traditionally, the spontaneous 
emission rate from the excited state Fq and the Rabi frequency 
n = El ■fJ'o/fi ifioll are specified in lieu of /Xq and El and we 
follow this convention here. We take Tq /2tt = AOAIHz in all 
that follows to model the organic dye terrylene in a hexade- 
cane Shpol'skii matrix at 1.7° K, a prototypical 2-level SMS 
system II19I1 . The following calculations assume various val- 
ues of A, ranging from 0.2 MHz to 200 MHz. Resolution 
down to 2 MHz is possible using a Fabry-Perot interferometer 
ll20ll . Theoretically, it should be possible to measure most of 
the reported quantities. 

A traditional measure of broadband photon statistics is 
Mandel's Q parameter IitIi. which is defined as the ratio of 
the second factorial cumulant of -/V(^ oo)(i) = N{t) to the 
first factorial cumulant (i.e the average) of N{t). Q (t) = 
[{N^{t)) - {N{t))^ - (Nit))] /{N{t)). We introduce a gen- 
eralization of this quantity appropriate to photon counting 
within a finite size frequency window 



it)) - (iV(^,A)(i))^ 



(^(.,A)(i)) 

(^(l,A)W)-(^(-.A)W)' 



(10) 



(^(-,A)W> 

and Qoo{^j->,t) = Q (t). Fig. [T] plots both the emis- 
sion lineshape (with finite resolution) I{uj — ujq) = 
limt^oo ^(-^(w.A) (^)) and Q A (t^, oo) for resonant excitation 
conditions (ljl = wq) and A = ro/4007r = 0.2 MHz. 
Two different values of the Rabi frequency are considered: 
il ^ To/ \/2 and = SFq / v^. The effect of frequency bin- 
ning is barely discernible in the lineshape when A is chosen 
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FIG. 1: Predicted normalized line shape (top panes) and (5a(w — 
iJott ~ oo) (bottom panes) for a 2-level dye with Fo /27r = 40MHz 
and two different Rabi frequencies Q, as indicated. A — 0.2 MHz. 
Insets truncate the y-axis to fully display the lineshape outside the 
vicinity of the high coherent peak at oj — ujq. 
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FIG. 2: QA{i^, t — oo) for different values of A and the same cases 
as in fig. [T] The bar graph format is used to emphasize the size of and 
location of the frequency bins, but is absent at the finest discretization 
for clarity. 



so small. Our results are essentially identical to the classical 
emission spectrum of Mollow lfl6ll . excepting the delta func- 
tion "coherent" iflollTill contribution at = ujl, which adopts 
a finite height after frequency binning. Plots for Qa oo) 
have not been reported previously, and at first sight our results 
appear surprising. The values selected for il in the chosen ex- 
amples both yield sizable negative values for the traditional 
broadband Q parameter (—3/4 and —0.11 for = Fo/V2 
and il = 5Fo/a/2 respectively), however (3a(w,oo) is seen 
to be positive over the entire frequency axis. The implication 
is that the antibunching phenomenon associated with Q < is 
due to correlations between photons of different frequencies. 
To make this point more explicitly, we plot Qa (w, oo) for dif- 
ferent choices of A in fig. |2] Qa{^^, oo) is seen to become 
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FIG. 3: Contour plot of the normalized factorial covariance function 
(see eq. [TTt . We consider the $1 = 5r/\/2 case of figs. [T]and[2]and 
have set A = 2 MHz. 

negative over portions of the frequency axis as A approaches 
the width of the peaks in the spectrum. Related behavior has 
been predicted for the intensity correlation function {g2) of 
photons originatingfrom a single well-resolved sideband in 
the Mollow triplet Isllgl [Toll . In that case, antibunching may 
be explained via the allowed sequence of photon emissions 
in the radiative cascade predicted by the dressed atom pic- 
ture ifioll . Interestingly, narrowband bunchinghas previously 
been attributed to properties of the detector Q, but we find 
the same effect in our observables that focus solely on photon 
emission. 

Eq. [8]is easily generalized to calculate correlations between 
photons at different frequencies. We define a normalized pho- 
ton covariance function as: 

CA{(^i,(^j,t) = (11) 
(iV(^..A)(t)iVK,A)(t)) - (iV(^..A)(t))(iV(c.,.A)(t)) ^ 

y^{iV(^.,A)(t))(iV(„,,A)(t)) 



A discretized version of this correlation function in the limit 
t ^ oo is plotted in fig. [3] where uj^j-^ have been chosen to 
follow LUi = ujQ + rA, where r is any integer. When oji = uij, 
CA(Wi,w.i,i) = (5A(t^i,i)- Otherwise, Ca(wj,Wj, t) sim- 
ply represents the covariance in photon number, normalized 
so as to give a finite result in the long time limit. Fig. [3] 
demonstrates that although QA=ro/40ir('^ii oo) > 0, the total 
Q parameter is dominated by negative contributions from pho- 
tons that are well separated in frequency; broadband measure- 
ment of Q contains important contributions from correlations 
spanning the entire spectrally active region of the transition. 
The positive inter-sideband peaks in fig. [3] reflect the corre- 
lated emission of photons from opposite sidebands. This is 
in qualitative agreement with the inter-sideband bunching ex- 
pected for a 2-level system excited far from resonance fst]. 
The phenomenon is attributable to the necessary paring of 
photons from the two sidebands in order to maintain total en- 
ergy conservation as photons of energy ftw^ are absorbed by 
the molecule. 

Our treatment of photon emission statistics is general and 
relies on no approximations beyond the RWA and Markov as- 
sumption for system dynamics. It is valid for arbitrary field 
strengths and does not assume particular physical regimes 
for the molecular system. Moreover, the present approach 
provides photon correlations between all possible frequency 
pairs, which enables calculation for any possible detector 
bandwidth and a quantitative demonstration of how seemingly 
inconsistent broadband versus narrowband statistics can arise 
from the same physical phenomena. This framework should 
prove valuable in the interpretation of future SMS experi- 
ments where moments higher than 1 will be measured and 
in understanding the molecular dynamics that such measure- 
ments probe. Several multi-state dye models are discussed in 
ref. 10] and will be treated in a future study ifTsll . 
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